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I. INTRODUCTION 



In the context of the study of QCD bound states via analytic methods a lot of interest 
has been devoted in the last ten years to the so-called Feynman-Schwinger formalism |]T|- 
M. The main feature of the formalism is that it allows to write the 4-point Green function 
(at least in quenched approximation) only in terms of a quantomechanical path integral over 
the quark trajectories times a functional depending on the average over the gauge fields of 
the Wilson loop defined by the quark paths. Moreover this functional can be expressed in 
terms of path derivatives of the averaged Wilson loop @J^]. Once we assume an analytic 
behaviour for the Wilson loop, which up to now can be only given as an external input more 
or less motivated by QCD but not completely derived from QCD, the advantages of such a 
formulation are apparent. It permits numerical calculations which under some conditions 
can provide results faster and cheaper than a traditional lattice calculation. Moreover it 
allows analytic estimates of physical interesting quantities. In Pj3|,f7|,|9yiOf it was possible in 



this way to obtain the complete heavy quark potential up to the order 1/m 2 for different 
Wilson loop assumptions and reproduce also the spin-dependent contributions of Eichten and 
Feinberg |Tl[] in the appropriate limits flIU |. We can say that the complete semirelativistic 



heavy quark-antiquark dynamics (at least in the form of the interaction potential) could be 
accessed only using this Feynman-Schwinger and path integral formalism. 

On the other side the derivation of the relativistic quark-antiquark interaction is a long 
standing and important problem (see j|] for some report papers). All the calculations of 
phenomenologically relevant quantities such as the masses and the form factors for the light 
hadrons rely on the understanding of relativistic quark dynamics [] . In the last years a lot 
of effort has gone into the development of light cone Hamiltonians on one side or Bethe- 
Salpeter-like and/or Schwinger-Dyson equations on the other. Some criticism has been 



1 In this paper we take into account only analytic models of the quark dynamics. Of course QCD 
lattice calculations are an alternative and complementary approach to the problem. 



made to the latter approach essentially related to the loss of gauge invariance |12[]. Indeed 
the manifestly gauge invariance of a physical state is a relevant concept when dealing with 
non-perturbative QCD dynamics. It is clear that the propagator of a coloured object could 
not be considered separately from the other coloured partners since it is connected by a 
string to it and this confinement dynamics dominates at large distances. Another way to 
put the thing is to say that non-perturbatively the background fields and their effect on 
the quark dynamics are important. Nevertheless we believe that when the gauge- invariance 
issue is properly addressed (i. e. the average on all the vacuum fields in the amplitude is 
correctly handled) the resulting effective interaction can be still treated in the framework of 
the Bethe-Salpeter equation and this supplies us with a formidable tool for the (numerical) 
evaluation of a huge number of physical quantities 0. 

One of the usually claimed limitations of the Bethe-Salpeter approach is that the con- 
fining part of the kernel is not known. In the literature it is widespreadly used a kernel 
made of a one-gluon ladder short range part plus a long range confining part suggested by 
a trivial relativistic generalization of the static linear potential. This amounts to consider a 
kernel depending only on the momentum transfer Q and with the form l/Q 4 . The Lorentz 
structure of the confining kernel is suggested to be a scalar again on the basis of the po- 
tential (which actually pertains to a complete different dynamical region, we point out) or 
simply phenomenologically treated as a vector which is chirally symmetric. However, all 
these assumptions run into great conceptual and concrete difficulties and it emerges that 
the kernel should be more complicate that a pure convolution type fli]. Another approach 



deals with Bethe-Salpeter and Schwinger-Dyson "coupled" equations with a kernel inspired 
by the lattice evaluation of the gluon propagator in a given gaugeQ |0J. 



2 For an example of the application of the Bethe-Salpeter equation to several phenomenological 
quantities see e.g. ]l3| . 

3 In this approach the main features of the phenomenology connected with the chiral symmetry can 
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The main motivation of this paper is to investigate the nature of the fully relativistic 
quark- ant iquark dynamics in the form of a Bethe-Salpeter kernel working with the Feynman- 
Schwinger representation of the quark- ant iquark gauge-invariant Green function []. The idea 
is to use this representation, that displays the complete dynamics factorized in the Wilson 
loop, in order to enforce the information we have on the Wilson loop behaviour directly on the 
Bethe-Salpeter kernel by means of a completely relativistic and non-perturbative procedure. 
This means that, starting with a form for the Wilson loop we are able to establish the 
leading Feynman graphs that make up the interaction kernel. Moreover if we use a Wilson 
loop behaviour containing the relevant part of the confining dynamics we will end up with 
the relevant part of the confining kernel. As it will become clear, the task is not simple in 
the case of quarks with spin. 

In practice a good part of the paper is devoted to the technical setting up of the formalism. 
As an application we derive the leading binding contribution to the Bethe-Salpeter kernel 
in QED (the one photon exchange graph). This supplies us with a technique and a definite 
language to apply in QCD. The up to now available assumptions on the behaviour of the 
Wilson loop average seem not to allow easy extensions. Therefore, we suggest the use of 
the Fock-Schwinger gauge in order to implement (as in the QCD sum rules approach) non- 
perturbative physics in the Wilson loop, leaving the structure of the Wilson loop average as 
close as possible to the QED one. We indicate the graphs relevant to the quark-ant iquark 
binding that make up the interaction kernel. 

The paper has the following structure. In section 2 and section 3 we derive the 2-point 



be qualitatively reproduced using a generic infrared enhanced gluon propagator. This is a further 
motivation of our believe that the characteristics of the light mesons can be well understood in a 
Bethe-Salpeter framework. 

4 Several attempts have been made also recently to obtain in the Feynman-Schwinger formalism 
a Bethe-Salpeter kernel for the QCD bound state H, but the problem is still open. 
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and 4-point Green functions in the Feynman-Schwinger formalism. In section 4 we apply 
the formalism to QED and in section 5 we discuss QCD and draw some conclusions. 



II. THE FEYNMAN SCHWINGER REPRESENTATION OF THE FERMION 

PROPAGATOR 

The aim of this section is to represent in terms of a quantomechanical path integral the 
fermion propagator S of a particle m in an external gauge field A. We assume A to be the 
non-Abelian gauge field associated with the gluon in QCD. Therefore where necessary we 
explicitly denote with the symbol P the path ordering prescription. In the Abelian case this 
prescription is obviously not needed. 

S is defined as Q 

S a/3 (x, y; A) = -(Tif) a (x)ipp(y))A , (2.1) 

where the brackets ( ) a stand for the average over the fermionic fields in the presence of the 
external source A. S satisfies the equations: 

(ip x -m)S(x,y;A) = 5 4 (x-y), (2.2) 
S(x,y;A){i % +m) = -6\x - y) , (2.3) 

where D% = d% — igA^(x) and D x = d x +igA^{x). 
If we define 

(z p x + m) A(x, y; A) = S(x, y; A) , (2.4) 

or alternatively 

5 All the formula here and in the following are given in the usual Minkowski metric (for our 
purposes we do not need to introduce the Euclidean metric which, of course, would be necessary 
in a more formal discussion). 
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A{x,y-A){-i p y +m) = S(x,y; A) , (2.5) 
then the function A satisfies the equation: 

(p x 2 + m 2 ) A(x, y; A) = -5\x - y) , (2.6) 
which can be written after some algebraic manipulations as 



-D 2 X -m 2 + -ga^F^x) ) A(x,y; A) = 5\x - y) , (2.7) 



1 

r 

with F^y = i [Dp, D u ] / g and a**" = i [7^, 7"] /2. In what follows it is useful to introduce the 
operator H(x, d x ) = (D 2 . + m 2 ) /2 — g a^ v ' F^ u (x) /A. Therefore, Eq. ( p.7| ) can be written as 

- 2 H(x, d x ) A(x, y; A) = 5\x - y) . (2.8) 

Following fl[ we consider the equation 

U A - H(x, d x )^j $(x, y; A;T)=i 8{T - T ) 5\x - y), (2.9) 

with boundary condition Q(x,y;A, T) = for T < T . The parameter T is usually called 
proper time. Eq. ( |2.9| ) is Schrodinger-like. The solution can be written as a path-integral 
over all the trajectories joining the point y at time T and the point x at time T (see e.g. 
0): 

$(x )2 /;A;T) = 0(T-T o ) Z(x,y;T; A) , 

rx=z(T) i / dt pz — H(z,p) 

Z{x,y;T;A) = VzVpVe J ?a 

Jy=z(T ) 

POO 

Integrating Eq. ( |2.9| ) in / dT and taking in account that Z(x,y;T ; A) = 5 4 (x — y), we 

J To 

obtain the solution of Eq. (|2.8| ) as 



1 f 00 

A(x,y;A) = -- dT $(x, y; A; T) 

2 JTo 



To 

T 

00 px=z{T) I 



fT 

/ dt pz — H(z,p) 

JTo 



. dT VzVpPe J To . (2.10) 

2 J T Jy=z(T ) 



Since the dependence on the momenta is Gaussian, the explicit integration on p is possible. 
Performing it we obtain 

r'l' jp 1 _j_ | 

if 00 , r =z(T) ~ { \ dt o gA"{z)z ll - 7 g<T^F lu ,{z) 

A{x,y;A) = — / dT / VzPe J t 2 4 

2 Jt Jy=z(T ) 

(2.11) 

Eq. (|2.11|) with Q2.4|) or (|2.5|) supplies us with a path integral representation of the fermion 
propagator in external field. We call this representation the Feynman-Schwinger represen- 
tation of the fermion propagator (some historical references are given in [|l~7|j ). 

III. THE 4-POINT GREEN FUNCTION AND THE WILSON LOOP 

Let us now consider a fermion-antifermion system. The corresponding 4-point Green 
function (see Fig. II) is given by 



1 r - i d A xC(ib ib A) _ 

G(x 1 , x 2 , yi , y 2 ) = —. / V^VAe J >(^) ^(vi) <% 2 ), (3.1) 



Af 

where M is a normalization factor and £ is the Lagrangian density of the gauge theory which 
we are considering (in our case QCD; in the following £ym will denote the Yang-Mills part 
of this Lagrangian density). Since the Lagrangian is quadratic in the fermion fields, it is 
possible to perform explicitly the integration over it. Neglecting 

i) fermion loops (quenched approximation), 

ii) annihilation graphs, 
we obtain f|[| 

lr if d 4 xC YM (A) 
G(x 1 ,X2,y 1 ,y 2 ) = jj= J VAe J iS(x 1 ,y 1 ; A) iS(y 2 , x 2 ; A) 

= (tS(x l ,y 1 ;A)iS(y 2 ,x 2 ;A)}. (3.2) 
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In order to deal with gauge invariant quantities, we will consider in place of the above gauge 
dependent Green function, the so-called gauge invariant Green function G inv obtained from 
the previous one by connecting the end points with the path-ordered operator 



ig / dzTA^z) 

U(y,x;T yx ) = Pe M* , (3.3) 

where the integration goes over an arbitrary path T yx connecting x with y. Within the 
approximations i) and ii) we have 

G inv (x 1 ,x 2 ,y 1 ,y 2 ) = (Tr iS(x 1 ,y 1 ; A)U(y 1 , y 2 ; P yiy2 ) iS(y 2 , x 2 ; A)U(x 2 , X\\ r^J) . (3.4) 

Writing now the fermion propagators in terms of the Feynman-Schwinger path integral 
representation given in the previous section, we obtain 

rTi m 2 _|_ z,2 

1 / ,,x poo r x 1 =z 1 (T 1 ) —1 / dti 

G hlv (x 1 ,x 2 ,y u y 2 ) = -(TTP(t]p^ +m) / dT x \ V Zl e Mo 2 

rT 2 m 2 + z 2 r 

r oo r y 2 =z2(T 2 ) -I \ dt 2 in & dz^AJz) 

x dT 2 Vz 2 e Mo 2 e y / r M ; 



T20 JX2=Z 2 (T2o) 

* f Tl dh 9 -a^F^( Zl ) 1 f T2 dt 2 9 -a$F^(z 2 ) <-( 2 ) \ 
x e Mo 4 e Jt 20 4 +m ) \ ; 

(3.5) 

where the upper-scripts ^ and ( 2 > refer to the first and second fermion line. T is the closed 
loop defined by the quark trajectories Z\(ti) and z 2 (t 2 ) running from y\ to X\ and from x 2 
to y 2 as t\ varies from T 10 to T\ and t 2 from T 20 to T 2 , and from the paths T yiy2 and T X2Xl 
(see Fig. Q). The quantity 

ig i dz^AJz) 

W(T;A) = TrP e Jr M \ (3.6) 

is known as the Wilson loop flHf|. 

Let us make some general statements. The variation with respect to the path of the path 



ordered operator U is given by (see for example [|I~9|,[2"0H ) : 



■I 



5U(y,x;T yx ) = igP<5y fl A fl (y)U(y,x;Ty X ) - Sx^A^x) U(y, x; Y yx ) 



1 z p 5z v — z^Sz^ ) 
ds F ltt ,(z(8))U(y } x;Ty X )>, (3.7) 



where we have assumed the path T xy to be parameterized by the proper time s in such a 
way that z(0) = x and z(l) = y. From it we have immediately: 

SU(y,x;T yx ) 



(is7i „ M -igP{F^(z)U(y,x;T yx )}, (3.8 
5U(y,x;T yx 



s igP\ \(y) U (y, x; T yx ) - ds z p —F pX (z) U (y, x\ T yx ) j> , 

5U{v 5 l'f yx) = ig p { -M*) u (y> *\ r ^) - 1] ds iP j^ F ^ z ) u (v> x > } 



(3.9) 



(3.10) 



where 5S pu (z) = dz p 5z v — dz v 5z p is the infinitesimal area. 

Let us now go back to the Wilson loop ( |3.6| ). As a consequence of Eq. (|3.8|) the insertion 
of a field strength tensor F^ v on a point z of the loop V in presence of the Wilson loop W 
can be written as 

s Ysl{t) = ~ tgF {w{Ti A) F ^ )] • (3 - n) 

This is known as the Mandelstam relation. Let us now assume that the string T X2X1 is a 
straight line. This is always possible since the string is arbitrary. We parameterize T X2X1 as 
z^{s) = Xi + s(x2 — Xi) M . From Eqs. (|3.9| ) and (|3.10|) we have 

SU(x 2 ,Xi,T xx ) . f 

j^jl = igP^A^(x 2 )U{x 2 ,x 1 ;T X2Xl ) 

- dss(x 2 - xi) p F pp (xi + s(x 2 - x 1 ))U(x 2 ,x 1 ;T X2X1 )\, 

(3.12) 

5 U{X2, x\] T X2X ) . ( 

fi^i =*</P| - A fl {x 1 )U{x 2 ,x 1 ;T X2Xl ) 

-J tfa(l-s) (ar 2 -a;i) p ^(a;i + s(x 2 - x 1 ))U(x 2 ,xi,T X2X1 )\. 

(3.13) 



Therefore we have 

(Tr A(x x , y l5 A) • • • f/(x 2 , x l5 r xaSl )) 

= (Tr d^ l A{x 1 ,y 1 ; A) ■ ■ ■ U{x 2 , z x ; r^J) 

+ (TrPA(xi,yi;A)---f/(a;2,Xi;r ;C2X1 )(-2 5(A /1 (xi))) - / ds (1 - s) (x 2 - x x ) p 

JO 

x (Tr P A(xi, yn A) • • ■ U(x 2 , x x \ T X2Xl )i g F^(x x + 5(0:2 - Xi))>, 
and finally (taking also in account ( |3.8| )) 

(Tr D^A(xi, y x - A) ■ ■ ■ U(x 2 , a*; T X2Xl )) = 

" Jo ^ (1 " S) (X2 ~ Xl) '6S„{ Xl + 8{x 2 - Xl ))) 

x (TrA{ Xl , yi ; A) ■ ■ ■ U(x 2 , xf, T X2X1 )). (3.14) 

In an analogous way we obtain 



(Tr - • • A(y 2 ,x 2 ; A) D U(x 2 , Xl ;T X2Xl )) 



(Tr • ■ ■ A(y 2 ,x 2 ; A)U(x 2 , Xl ; T X2Xl )) 

lo ' " ~ 1Jt> SSp^xx + s(x 2 - xi)) 
Therefore Eq. (|3.5|) can be written as 



x I d x 2 ~ / dss(x 2 -x 1 ) p — — ■ — — )■ (3.15) 



x 2 ,yi,y 2 ) 

4 ~ ^ Jo v ^ " yv ~ z ~ iy ^^(Xi + s(x 2 - x x )) 

/■oo /•xi=zi(Ti) — Z / tXti 

x / dTt Vzxe j t w 2 

•Mo Jyi=zi(T w ) 



' <! >- -' Y Io dS (1 - S) {X2 - Xir 5S^ + 6 s(x,-xi)) + m 



(i) 



X 



rT 2 m 2 _|_ _j2 
00 ry 2 =z 2 (T 2 ) ~i / at 2 » 

dT 2 / D2 2 e 7t 2Q 2 

T20 ^2:2=22(720) 

2 



x« - T io 4 5^(zi) e Mo 4 ^^W^r.A)) 

(2) 

i?„+ir />'fa-».)' M , (J1+ < . (j ,_ Jl)) +"■) • 
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All the dynamical information are contained in the Wilson loop average (W(T,A)) and in 
its functional derivatives. The analogous happens in potential theory where it is possible to 
express the potential up to order 1 / m 2 only in terms of the Wilson loop functional derivatives 
0J^{TI| . If we were able to know exactly the Wilson loop average over the gauge fields, then 
we could express the 4-point quenched Green function as a pure quantomechanical path 
integral (which is very convenient also for numerical applications see for example ||). This 



would realize the Migdal program of [|2(J . Of course the difficult point is to give an evaluation 
of the Wilson loop. In the next section we will discuss the QED case, for which the Wilson 
loop average is exactly known in analytic closed form. In particular in order to see how Eq. 
( p,16| ) works we will show how to recover the binding interaction kernel in terms of Feynman 
graphs. 

IV. AN EXACT APPLICATION: QED 

Since in QED the Yang-Mills Lagrangian is quadratic in the fields, the Wilson loop 
average can be evaluated exactly. The result is 

— — i dx^ i dy u D uu (x — y) 
(W(T,A)) = e 2 It Jy ^ y \ (4.1) 

were D^ v (x — y) = (TA^x) A v {y)) and g has now to be interpreted as the electron electro- 
magnetic charge. In the quenched approximation D^ v is nothing else than the photon free 
propagator. Let us consider Eq. (|4.1|) only up to order g 2 . Therefore our expression for the 



Wilson loop in QED will be 

2 

(W(T,A)) = 1 -jfcdx" f T dv v D^{x-y). (4.2) 



Limiting ourselves to the Wilson loop expression (|4.2| ), the next step will be the evaluation 
of the area derivatives which appear in (|3.16| ). A simple calculation leads to 

= -9 2 j/y v (d,D av (z -y)- d a D Pv (z - y)) . (4.3) 
11 



As a particular case we have 

I ds * - s ^ I* - ^ rn^^^s) = ~ g2 i m ^ - y) 

g 2 £ dy v jT 1 ds d, [x p D pv [x + x 2 - y))| x=s(xi _ X2) , (4.4) 

~ g2 i dyU J* ds d " {xPDpAx +Xi ~ ^u**-** ) • ( 4 - 5 ) 



Since the contributions in the second line of Eqs. ( f4.4| ) and (|4.5| ) are exactly canceled by 
the action of the derivatives $ and (/) on the string T X2Xl , we will assume in the following 
that these derivatives do not act on the endpoint strings as well as we will take into account 
only the contribution of the first lines in Eqs. ( J4.4Q and ( |4.5| ). This will simplify the display 



of the results. Putting Eqs. (fO|)-(|3~5D in i^J^> we obtain 

G inv (x 1 ,x 2 ,y 1 ,y 2 ) = 
I(i^ i+m )W...(-i^ 2+m )^) 



-- f Tl dh I dyV, Qj^D^zi -y)-- A Ux» F dt 2 [Y, $J 2) D, u (x - z 2 ) 



/T10 JT 4 JF JT 20 

~T£ C dU Cdt2[YjJ 1) [YjJ 2) D lll/ (z l -z 2 ) 

IO JTio JT20 

~h jC ^ It! ^ [T " ^ ](2)[Y ' ^'f )D ^ Z2 - + m)(2) 

2 ( IT 

~ • • • [ - *7 M(1) jf difivfa - y) - 47 M(1) jH dfelV, ?U (2) £U*i - ^2) 

jC rfti[7 '' ^ ](1)jD ^ (a:i - ^ + m)(2) 



^(z ^ + m) « • • • I - i7"< 2 > j£ dxPD^x - x 2 ) 



- V 2 ) r d tl [Y, $J 1)D ^ - x 2 ) - \l v{2) r dt 2 [Y, $ z f ] D» v {z 2 - x 2 )\ 

4 4 JT20 J 

+ j---7 M( V (2) ^(^i-^), (4.6) 
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where the dots indicate the path integrals and the kinematic factors given in ( |3.16| ), precisely 



T i m 2 + z\ 

!El=*l(Ti) -% I dti 

dTi I Vz x e Mo 2 

Tio Jyi=zi(T 10 ) 



ry2=z2(T 2 ) -I I dt 2 - 

x / dT 2 / Vz 2 e Mo 2 

M 20 JX2=Z 2 (T20) 



T 2 ^ m 2 _|_ £2 



(4.7) 



In order to recover Feynman diagrams from Eq. ( f4.6|) , we have to throw away the 
endpoints string contributions. Because we have already taken into account the action of 
the derivatives on these strings by neglecting the second line contributions in Eqs. ( }4.4|) and 
( |1.5| ), this can be done without any problem. Of course in this way we lose gauge invariance, 
but Feynman graphs, as usually known, are not gauge invariant quantities. Furthermore 
in QED the manifest gauge-invariance of the two particle state is not a relevant concept. 
Moreover we will neglect self-energy contributions (which for usual gauges do not contribute 
to the binding), following the replacement scheme: 



i dx^ i dy v 2 / dz£ / dz\ 
Jr Jr Jy 1 Jx2 

rv2 

dy v ^ dz v 2) 

I JX2 

dx^ r ' dzit. 

i J Vi 



Let us consider now the kinematic factors (fL7|). We define 



2 poo rx=z(T) — I 

A(x-y) = A(x,y;0) = — / dT Vz e Jt o 

2 Jt Jy=z(T ) 



dt 



z + m 



From the definition we have (see also Eq. Q2.71 )) 



A(p) = J d*ze ipz A(z) 



1 



pi _ m 2 _j_ j e ' 



and 



S(p) =i d 4 z (i$ z + m)e ipz A{z) 



•j) — m + ie 



Some properties of A are (see Appendix): 



(4.9) 



(4.10) 
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T z 2 + m 2 

r oo px=z(T) —i l at 

i) I dT Vz e Jt o 2 

J T Jy=z(T ) 

-4jd^A(x-0f(0^-y), 

f T z 2 + m 2 
poo rx=z(T) —i / at 

ii) dT Vz e j t 2 

Jt Jy=z(T ) 



T 



dtf(z{t)) 



(4.11) 



To 



dtf(z(t),z(t)) 



-4 J d'i J A J ^ e -^A(x-Of(^,p) A( V -y). (4.12) 

Since all the calculations are more transparent in the momentum space, it is convenient 
to deal with the Fourier transform of the Green function G, defined to be (see Fig. |]): 

G( P ' 1 ,P2,Pi,p' 2 ) = Jd 4 Xl Jd 4 x 2 Jd 4 Vl Jd%e-^e^y 2 e^e-^ 2 G(x 1 ,x 2 ,y u y 2 ). 

(4.13) 

The first term of the right-hand side of Eq. ( [4.6| ) is therefore nothing else than the 2- 
particle free propagator. At the order g 2 we have to evaluate the next four terms (we are 
considering only interaction terms). These four terms can be factorized out in the product 
of two contributions on the first fermion line and two contributions on the second fermion 



line. Taking into account only the first fermion line and using Eq. (^4.12|) we obtain from the 
first contribution a term like 

d A p 



-2 J d A Xl J d^e^e-^^^+m)^ J d^ J d\ J 



(2^) 4 

) &(7i-Vi) 



which after a straightforward calculation ends up to be 



- *f®DM -Pi--- -)A(pi) + ISip'^Sipi)] (1) D,M - Pl - ■ ■ •) 



(4.14) 



where D^ u is the Fourier transform of the photon propagator. The second contribution on 
the first line is 

d 4 Xl [ d% e^ e -^A(z a - yi)Y {1) D^( Xl -■••)= -f {1) I>M - Pi - • ■ -)A(pi). 



(4.15) 
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Combining together ( |4. 14| ) and ( 4.15 ) we obtain the usual fermion-photon vertex on the first 
fermion line: 

5(Kb^(Pi)] (1) d^(p[ -vx -•••)• 

Handling in the same way with the second fermion line we finally have 



G(p[,p 2 ,Pi,P 2 ) = (2tt) 4 5 4 (K - Pl )S(p' 1 )W (2ny6\p' 2 -p 2 )S(p 2 



I (2) 



+ (2tt) 4 5 4 (K +p 2 -pi- ft ) [5(pi) ^7 M S(P 
X-D^bi - Pi)- 



(i) 



s(P2)ig r r v S(p 2 [ 



(i) 



(4.16) 



This is the 4-point free Green function plus the one-photon exchange graph. We notice that 
in order to obtain the Bethe-Salpeter equation for QED all the higher order corrections to 
the formula (|4.2j ) have to be taken in account. This is a little bit cumbersome (and is beyond 
the purposes of this discussion) but can be done in a systematic way with the methods given 
in || for the scalar QED case. Here we want to point out that starting from the Wilson loop 
behaviour given in Eq. ( |4.1| ), and restricting ourselves to the g 2 contributions, we were able 
to reconstruct with this technique the ladder kernel. In a similar way, starting with a given 
confining behaviour of the QCD Wilson loop we should be able to reconstruct the relevant 
contribution to the quark-antiquark Bethe-Salpeter kernel. 



V. QCD CONFINING MODELS AND CONCLUSIVE REMARKS 

In principle the same technique can be used in order to extract an interaction kernel from 
Eq. (|3.16|) in QCD. The problems arise from the fact that we do not know the exact analytic 
expression for the Wilson loop average in QCD and therefore we do not have the equivalent 
of equation ([4.1] ). We have to resort to some model approximations and from this respect it 
is really determinant, in order to recover a meaningful result, to deal with a gauge- invariant 
quantity as the Wilson loop is. However the up to now available approximative expressions 
for the Wilson loop average in QCD seem to be too rough for this purpose. In fact in order 
to recover a reasonable interaction kernel from the Feynman-Schwinger representation of 

15 



the 4-point Green function, the dynamics of the interaction, contained in the Wilson loop, 
cannot be anything, but it has to fit with the kinematics of the two quarks, represented 
in Eq. ( |3.16| ) by the external and spin projectors and the kinetic energy terms inside the 
path integral. For example, this matching condition manifests itself clearly in QED by the 
cancelation of ( |4.15| ) when combined with ( 4.14 ). This is the difficulty when particles with 



spin are taken into account. It is a consequence of having expressed the field insertions in 
the Wilson loop in terms of functional derivatives on the Wilson loop contour. In this way 
not only the gluodynamics, but also the coupling of quarks and gluons (in the perturbative 
and non-perturbative regime) is depending from the Wilson loop behaviour. 

The authors of ref. || assume the Wilson loop average to be governed by the Wilson area 
law, i. e. 

(W(T,A)) = e - aS ^ T \ (5.1) 

where S , min (r) is the minimal area enclosed by the closed curve T and a m 0.2Gev 2 is the 
string tension. The matching between this dynamical assumption and the kinematics of the 
quarks is successful only in the so-called second order formalism. In other words, only if part 
of the kinematics is not taken into account. 

Also more sophisticated assumptions for the Wilson loop average are difficult to handle 
in the Feynman-Schwinger framework without making semirelativistic approximations. In 
the stochastic vacuum model (SVM) [||,^l|] one assumes that 



,2 



- 9 - [ dS^(u) [ dS^WiF^up^v^Fpxivpiv,^^)) 
(W{T,A)) = e 2 J SF) Mr) , (5.2) 

where the curve T uv connecting the points u and v is arbitrary and the integration is per- 
formed over a surface S(T) enclosed by the curve T. For what concerns the present discussion 
we can neglect color indices (the bracket ( ) is an identity matrix in colour space). In the 
usual straight-line parameterization of the surface (see for example |10|| ) we introduce points 
belonging to different fermion lines evaluated at the same proper time which seem not to be 
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treatable in the previous formalism. More convenient appears to parameterize the surface 
in triangle having a fixed vertex and two vertices running on the curve T. This corresponds 
essentially to choose the gauge fields in the so-called Fock-Schwinger gauge. We will discuss 
this case in more detail in the following. 

Let us assume that the fields satisfy the gauge condition: 

(x - x l YA lx {x) = 0. (5.3) 

As a consequence we can express Ap in terms of the field strength tensor Fp V f22|| : 

AJx) = / daa{x — xi) p FpJx\ + a{x — x\)). 
Jo 



This gauge is a very natural tool in the sum rules method [p^] and a great deal of the existing 
information on non-perturbative QCD dynamics can be recovered working in it. 

Expanding the Wilson loop average in cumulants and taking only the second order ones 
|2T| , we obtain an expression formally equivalent to Eq. ( [4. 1| ) where Dp V is now no more a 
local quantity (the gauge breaks in fact the translational invariance) and is defined to be 

Dp U (x,y) = daa(x- x^ J d(3(3(y- Xi)° (F m (xi + a(x - x^F^xi + f3(y - xi))). 

(5.4) 

Perturbative and non-perturbative contributions are contained in (F P pF all ). We focus our 
attention on the non-perturbative ones which are of the type: 

{Fpp{u)F au {v)) = (g pa gp U - g pi ygpa)f(a 2 (u - v) 2 ), (5.5) 

where 1/a w 0.2 -4- 0.3 fm is the correlation length which defines the confining energy (dis- 
tance) regions. Notice that in the limiting case a — >• the form factor / coincides with the 
gluon condensate. In this way we are considering the Wilson loop behaviour given by the 
stochastic vacuum model. The model PJ24]| is based on the idea that the low frequency con- 



tributions in the functional integral can be described by a simple stochastic process with a 
converging cluster expansion. Assuming the existence of a finite correlation length 1/a linear 
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confinement is obtained. The simplest formulation is characterized by a Gaussian measure 
specified by the correlator given in Eq. ( |5.5| ) just determined by two scales: the strength of 
the correlator (the gluon condensate) and the correlation length. This behaviour of the cor- 
relator (i. e. the exponentially falling off, like a gaussian, of the function / as \u — v\ ^> 1/a 
in Euclidean space) has been directly confirmed by lattice calculations [25]. The Wilson loop 
behaviour given by Eqs. (|5.4|) , (|5.5|) has been successfully in applications to the study 

of soft high energy scattering |M] as well as to the heavy quark potential. In particular the 



SVM potential reproduces exactly static |18| , spin-dependent |TT[ and velocity-dependent 
potentials @J^] in the appropriate limit jlQfl . 

Using, now, Eqs. ( }4.3|) - ([4.5|) we obtain Eq. (|4.6|) but with the above definition of D^. 
We observe that because of the gauge condition (|5.3|), either the second line of Eq. ( |4.4| ) 
and Eq. ( |4.5| ) vanishes either the action of the derivatives on the string T X2X1 is zero, being 
U(x2, Xi] T X2Xl ) = 1. If we assume that also in this case the endpoint string U (2/1 , 2/2; Fyiy 2 ) is 
not relevant for the bound state [] the Feynman graphs contributing to the interaction kernel 
are given in Fig. ^. With the box we indicate the translational non invariant propagator 
given in Eq. ( |5.4| ). In particular, due to the losing of translational invariance, the former 
(in QED) "self-energy" graphs, which now can be interpreted as the the interaction of the 
single quarks with the background vacuum fields, contribute to the binding and can not be 
neglected anymore. This last point emerges very clearly in the one body limit. If the mass of 
the particle moving on the first fermion line goes to infinity, then the exchange graph of Fig. 
|3] does not contribute at all to the static limit which is entirely described by the interaction 



6 In this way we lose gauge invariance. Nevertheless, this assumption seems to be justified by the 
existence of a finite correlation length in the correlator dynamics ([5.5|). Therefore, in the limit 
(T — To) — ► 00 the contribution of the string should be negligible. Anyway we stress that all 
the dynamics approximations have been made on gauge invariant quantities. Finally, in the case 
Ui = 2/2 the contribution of the string actually vanishes. 
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of the second quark with the vacuum background fields [57] . Moreover in the one body-limit 
we have shown |27J that in this kind of graphs the confining dynamics contained in the 
correlator combines itself with the quark propagator in such a way that the pole mass turns 
out to be shifted by an amount a. Eventually, in the limit a — > oo the quark propagator 
has no pole mass. This means that by cutting the Feynman diagram you could not produce 
a free quark at least for some values of the parameters, which seems to be in the line of 
the results obtained by the groups working with Bethe-Salpeter and Feynman-Schwinger 



equations with phenomenological kernels |15[]. Work is currently going on to further clarify 
this point. 

Finally, we observe that the Lorentz structure of the obtained kernel is not simply un- 
derstandable in terms of a vector or scalar exchange. 
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APPENDIX: 



In this appendix we will prove Eq. ( [4.12 ). Suppose to be e the proper time interval in 
the path integral definition: 



27V-2 

Vz 



( — ) d%2)...d±z{N-l). (Al) 

The path integral measure satisfies the following properties (2 < n, n + 1 < N — 1): 
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x r r z ( n ) r x 

Vz = / d 4 z(n) / Vz / Vz, (A2) 

y J Jy Jz(n) 

_ D , = /^(„)/d«z( B + l)(— ) / ^/ ( „ +1) ^, (A3) 
and in the path integral we have 

/w),iM) S /f ' (t)+ f (t +V (t+e) -* (t) V (A4) 



Therefore the left-hand side of Eq. ( |4.12j ) can be written as 

f T z 2 + m 2 

r oo r x=z(T) —i / Ctt rT 

dT Vze Jt 2 / dt f{z{t), z{t)) 

JTo J v=z (Tn) JTo 



T Jy=z(T ) J T 

. /•* z 2 + m 2 

r r / ^ \ 2 r°° r T rv=z(t) —i I dr 

/>£ d ± v (J_) dT dt Vze Jto 2 

J J \2meJ J T J T Jy=z(T ) 

. f T , z 2 + m 2 rt+e , z 2 + m 2 

r x=z(T) —i I dr — i / dr 

x / Vz e J t+£ 2 e Jt 2 

x/ Mt)+z(t + e) z (t + £ )-^(t) \ (A5) 



In the e — > + limit we have 

2 -i J dT — J <t)+<t + *) z(t + e)-z(t) 



t+e j ^ _|_ m ^ 
C 

/ 1 \ 2 ~* v " v " ' ~i — , fz(t) + z(t + e) z(t + e)-z(t)' 



2 /, 

/ V 2 

; (,(t + e)-,(t)) 2 



V2W J \ 2 e 

, f z(t) + z(t + e) , d \ — 

J I 2 ' dz(t + e)J \2me) 6 



I z(t) + z(t + e) d \ r d A p 
f [ 2 ,l dz(t + e) J 



-ip(z(t+e)-z(t)) 



dz(t + e)J J (2tt) 4 
d*p ip(z(t+e) _ z(t)) ( z(t) + z(t + e) 



(2< 



/*oo /*J /*oo /*oo 

Putting now (|A6| ) in Eq. (|AI|) and taking in account that / dT / dt = dt dT we 

-'To «^To J To •/ t 

obtain Eq. ( |4.12j ). As a particular case of Eq. ( |4.12| ), if / does not depend on z, Eq. ( [4.1 1| ) 



follows immediately. 
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FIGURES 




FIG. 1. The 4-point Green function. 
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FIG. 3. Non-perturbative contributions to the bound state kernel in QCD. 
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